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Abstract 

The radiative meson decays V — > -P7 and P — > 77 are analyzed using the 
quark triangle diagram. Experimental data yield well determined estimates 
of the universal quark- antiquark- meson couplings gvqq 1 and gp q q' for the light 
meson sector. Also predictions for the ratios of neutral to charged heavy me- 
son decay coupling constants are given and await experimental confirmation. 
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I. INTRODUCTION 



In earlier work |l| we used a supermultiplet theory uniting the vector and pseudoscalar 
mesons to attempt to obtain a universal three-point coupling constant. The relativistic 
multispinor fields of the supermultiplet theory described pointlike mesons with correct spin, 
parity, flavor and color degrees of freedom without necessarily invoking the notion of con- 
stituent quarks. Despite the apparent contradiction with the modern understanding of the 
quark nature of mesons it has been shown M that such a field is dynamically equivalent 
to a system of two quarks moving at equal velocity and on-shell. This relatively simple 
scheme does nonetheless compare favorably with experimental results for strong vector and 
pseudoscalar interactions in the light and heavy meson sector. When we include the radia- 
tive decays by incorporating vector meson dominance with the scheme we once again found 
reasonable agreement, but there were some unexpected and significant deviations from the 
theory, particularly in the K* — > K'y decays. 

The supermultiplet theory and exact SU(3) predict the coupling ratio IgK^K^/ 9k*+k+^\ 
to equal 2. But experimental measurement currently estimates the ratio as 1.51 ± 0.13, a 
substantial difference. One possible reason why the supermultiplet scheme did not comply 
with the experimental measure is that the exact form of the vector meson dominance is not 
known in the q 2 — > limit; it is only accurately known for q 2 = rriy from V — > decays. 
Hence there is some uncertainty in extrapolating vector meson dominance to the off-shell 
case. However, the difference in theory and experiment is so large that it is unlikely this 
is the only contributing factor. Thus we also implemented some symmetry breaking in the 
supermultiplet scheme, but the K* radiative decays seemed impervious to our attempts at 
matching theory with experiment as large discrepancies remained. 

To understand these deviations further, we sought a method which easily allowed for off- 
shell propagation of the quarks so that we could evaluate the magnitude of this necessary 
correction. A convenient and apparently successful method for doing this was by use of a 
quark triangle diagram; which so far has given accurate predictions for 7r° — > 77 decay widths 
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and pion and kaon charge radii ||[|]. A form which used chiral and isospin symmetry 
has also been successfully applied to the K* radiative decay problem and some radiative 
decays in the light meson sector || . The resulting loop integral accounted for the difference 
in quark masses and hence propagators in the loop, and correspondence between theory and 
experiment was achieved. 

Nonetheless, a crucial assumption of the quark triangle diagram is that the meson-quark- 
antiquark vertex has the form gpgq' 75 for the pseudoscalar meson and gvqq'l' 1 f° r t ne vector 
meson. If we are to confidently use the quark triangle method we wish to test the appro- 
priateness of this assumption. We do this by extracting the couplings from experimental 
measurements and examine the extent to which they carry the spin and flavor symmetries. 
With this in mind, we firstly formulate the limit free form of the integral and derive a limit- 
ing case which uses chiral symmetry. This enables us to compare our result with others. We 
also determine a heavy quark expansion. Following this, the scheme is applied to P — > 77 
decays which yield estimates of gp q q> and then to V — > P7 decays to obtain the product 
Qvqq'QPqq' f° r different quark flavors. 

The results indicate that the meson-quark-antiquark couplings in the VVP sector de- 
termined from different channels which involve common constituent quarks are remarkably 
uniform, suggesting that the effective vertex in the quark triangle diagram is valid. The 
data also demonstrates that the triangle method should be highly predictive due to the 
stability of the couplings. Finally we use the method in the heavy meson sector, to predict 
coupling ratios of the form gyopo^/ g v + P + 1 where the only free parameters required are the 
constituent quark masses. Our result for D* decays fall within other theoretical estimates, 
while that of B* is sensitive to b quark mass. 

II. THE QUARK TRIANGLE 

In the quark triangle formulation of Figure [I], the decay from vector meson to pseu- 
doscalar meson and photon state is mediated by a quark loop with flavors of constituent 
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mass m and m. (The choice of constituent mass rather than current mass is supported 
The quark triangle diagrams correspond to the Feynman amplitude for the decay, 



A(V -> Pj) = -N c gv q q>gp q qieQ^e v 



fo \4 Tr ^^ 7 /, ^ — 52 -^jTH ] + (m ^ m,Q ^ Q), (1) 

(Z7r) 4 p — m p+pA*—f — m p— jc — m 

where K^{e v ) is the vector meson (photon) polarization vector, M*(M) is the vector (pseu- 
doscalar) four-momenta, gvqq'(gpqq') is the vector-quark (pseudoscalar-quark) coupling con- 
stant and eQ is the electric charge of the quark of mass m in the loop. 

The Feynman loop integral involved in the amplitude (0) may be solved with standard 
techniques. We maintain the notation of Bramon and Scadron || and call the integral J, a 
dimensionless quantity after multiplication by m. Subsequently, 

t f 1 t f 1- " r m + (m — m)u , . 

J = ml du I dv — -—^ —rr, , 9 7z r (2) 

mr + ywr — m z )u — M* z uv — M z u[l — u — v) 



o Jo 

m 2 n 



( du{5 + l/u){\n\f\-\n\j\) (3) 
Jo 



M* 2 - M 2 

■[Jl-Jx + JZ-J, 



2 

m 



M* 2 - M 2 



where 



m — m 
o = 



m 

j = [ m 2 - (m 2 - m 2 + M 2 )u + M 2 u 2 ]/m 2 

J\ — — I duhi\j\/u 
Jo 

J 2 = — 5 / du In Ijl, 
Jo 

and j*(J*, J|) corresponds to j(Ji,J2) with M M*, respectively. Now M*(M) corre- 
sponds to the vector (pseudoscalar) mass. 



A. Determination of Ji 



In attempting to find an expression for J 1; we re- write the argument of the natural log 
in a form similar to that of the dilogarithm. To do this, we factorize j as 
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j = 1 + (S 2 + 25 - (M/mf)u + M 2 u 2 /m 2 = (1 - v lU )(l - v 2 u), 

where 

Vl>2 = -{5 2 + 25- (M/m) 2 + [(5 2 + 25 - (M/m) 2 ) 2 - (2M/m) 2 } 1/2 }/2 (4) 
= {m 2 -m 2 + M 2 ± A 1/2 (m 2 , m 2 , M 2 )}/2m 2 . (5) 

and 

A(m 2 ,m 2 ,M 2 ) = [M 2 - (m + m) 2 }[M 2 - (m-m) 2 ]. (6) 

We obtain a similar expression for u* 2 upon substitution of M by M*. The factorization we 
have performed does not necessarily lead to real v* k or and we must consider the case for 
both real and complex arguments. 



1. Real vl or 

We first consider the simplest case, that when either vl or Vk is real. For the moment 
we simply deal with real Vk and extend our findings to real by substitution of M by M*. 
The factor Vk is only real when A(m 2 , m 2 , M 2 ) > which from (^) implies 

M > m + m ot M < \m — rh\. (7) 

When we are assured of real Vk the solution of J\ is related to the standard dilogarithm 
function, 

j, = _ /' du ± hli^A ^ £ U2(Vb o) = ( Ll2(t ") for 5 1 (8) 

° fc=1 M fc=1 [ ELiLi 2 (^) + ^ln|^| for v k > 1, 

Complex vl or v k 

The ffc are complex if A(m 2 , rh 2 , M 2 ) < and we need an appropriate method for handling 
this situation. Fortunately the dilogarithm of a complex argument does exist so we may 
proceed. 



We express J\ as 



k=i u 



Ji 



k=l 



1 fp ln(l — 2xcos0fc + x 2 ) 



x 



Jo 



dx + i arctan 



V sin fc 
1 - y cos fc 



1/ 



Jo 



ln(l — 2a; cos + a; 2 ) 



x 



= 2Li 2 (p, 0), 



(9) 



where 



p = M/m, = 0x = —02 and cos0 = 



m 2 — m 2 + M 2 
2Mm 



and J* corresponds to J x with M <-> M* 



B. Determination of J2 



(10) 



J 2 is a simpler integral to evaluate as it does not contain the \ ju dependence. Recall, 

J 2 = -5 f 1 du lnfl - (1 - m 2 /m 2 + M 2 /m 2 )w + M 2 u 2 /m 2 ], 
Jo 



which, by standard techniques, reduces to 



J2 



m — m 
m M 2 



.. m 2 -m 2 + M 2 ) In - - A 1/2 (m 2 ,m 2 , M 2 )arctanh ( 



m 



and is valid for all m, m, M. We obtain a similar expression when we substitute M — > M*. 

It is also useful to express J 2 in terms of real v k - Following a similar method to that 
used for deriving Ji, we find 
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j 2 = Sj2(l-l/v k )\n\l-v k \ for A(m 2 , m 2 , M 2 ) > 0. 
fc=i 



(11) 
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III. COMPARISON WITH CO VARIANT AMPLITUDE 



Our final form for the loop integral is: 



J = 



m 



\J{-Ji + J* 2 -h]. 



(12) 



M* 2 - M 2 1 

Here we have not considered the imaginary part in J 1 which is irrelevant to the decay process, 



ELi Li 2 K 0) if A(m 2 , m 2 , M 2 ) > 

where v 1>2 = [m 2 - m 2 + M 2 ± A 1/2 (m 2 , m 2 , M 2 )]/2m 2 
2Li 2 (p,0) if A(m 2 ,m 2 ,M 2 ) < 

where p = M/m, cos0 = (m 2 — fn 2 + M 2 )/2Mm 



(13) 



and 



m — m 



> , r -., m \i/2/ 2 -2 .^n , (\ 1/2 {m 2 ,m 2 ,M 2 ) 

m — m + M ) In A ' (m , m , M arctanh — - — ^ 

; m V ; \ m 2 + m 2 — M 2 



m M 2 

Thus, our Feynman amplitude for the decay is 



(14) 



where J — > J when m <-> m (from the momentum crossed Feynman diagram) and (Q is the 
charge of the quark with mass m. We compare this with the general covariant amplitude 
for the process V — > P 7 

A(V - P 7 ) = ig^e^P 01 ^^ 
so that our quark triangle approach resolves the 5V P7 covariant coupling constant as 



g VPl = N c eg Vq q,g Pq qi[QJ/m + QJ/m]/47r 2 . 



(15) 



IV. P -> 77 IN THE QUARK TRIANGLE SCHEME 



We are interested in understanding the behaviour and obtaining actual values for the 
coupling constants gp qq i and gy qq i in the light quark sector so that we can use appropriate 

7 



estimates for these couplings in the heavy quark sector. To this end we can use the well 
documented decay data for V — > P7 in the light vector meson sector, as well as the decays 
P — > 77. These latter processes are particularly useful as they only involve the coupling 
gpqqi, and not the product gvqq'QPqq' as do the first case. Consequently, we must derive 
the amplitude for the decay of a pseudoscalar meson into two photons. This we may do 
by following a similar derivation as above, but it is much simpler to make the following 
substitutions in the V — > P7 amplitude (|i~4f) : 



M* — ► M, M — ► 0, g Vq q' -> gpqq>, and g Pq q> -> eQ, 

and since all pseudoscalar mesons involved in P — > 77 decays must be quark flavor singlets, 
m = fh, Q = Q, J = J. Subsequently ( |15D is reduced to 



g Pll = 2N c gpqq>e 2 [Q 2 J/m}/4: 



71 



(16) 



where 



J = i 



f2ELiLi 2 (^,0) if M>2m 

where v 1<2 = M[M/m ± (M 2 /m 2 - 4) 1 / 2 ]/2m 
2^Li 2 (p,0) ifO<M<2m 

where p = M/m, cos<ft = M/2m 



V. CHIRAL LIMIT 

The chiral symmetry limit is useful in the light meson sector, and gives our work di- 
rect comparison with that of Bramon and Scadron ||. The limit corresponds to a small 
pseudoscalar mass when compared to the vector mass, that is M* 2 » M 2 . Such a limit is 
entirely appropriate for the study of the radiative decays of K* mesons, and it is re-assuring 
to know that our J reduces to the Jm of in the chiral limit. 

The chiral limit, corresponding to M — > in J, enables us to use the real form for J\ in 
(I) as M < I m — fh\. 
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and 



v k {M -> 0) 



for k = 1 

-(5 2 + 25) for jfe = 2, 



from (01) so that 



J x = U 2 (v k , 0) = Li 2 (0, 0) + Li 2 (-5 2 - 25, 0) = Li 2 (-5 2 - 25, 0), 
fe=i 



J 2 = -5^(1-1/^,) hi |1-^.| 



k=l 



-5-5[l + l/(5 2 + 25)]ln|l + 25 + 5 2 | = -5 - 2 ^ + ^ In |1 + 5\. 



5 + 2 



Thus J in the chiral limit, which we denote as J CL , becomes 

,2 r 2 



J 1 



CL 



m 



M 



*2 



5+^[Li 2 (^0)-5(l-lK)ln|l-^|]... 



fc=i 



Ti 2 (-5 2 - 25,0) + ^^ln|l + 5| 



(17) 



where we have assumed v* k is real. This form may be simplified even further near the isospin 
symmetry limit whereby m ~ fh. In this instance we ignore 5 terms of order 2 and higher. 



Li 2 (-5 2 - 25,0) 
2(1 + 5) 2 



2 + 5 



ln|l + 5| 



Li 2 (-25,0), 

(5 + I) i n |i + 2( 5|, 



y l,2 



M 



{25 - (M*/m) 2 + [(25 + (M*/m) 2 ) 2 - {2M* /m) 2 } 1 ' 2 } /2 

1/2 



*2 



2m 2 



-5 + 



M 



*2 



2m? 



5 



M 



*2 



(18) 



Thus we find J CL incorporating isospin symmetry between quark flavors reduces to 



J' 



CL 



m 2 { 2 

\ 6 + EtM**, 0) - 5(1 - IK) - vl\] 



-Li 2 (-25,0) + (5 + -)ln|l + 25| 



with vl defined in Eq. flIB|). This form is very similar to the Jm of Bramon and Scadron 
||. There is a subtle difference in their use of the dilogarithm function Li 2 (z) versus our 
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function Li2(z, 0) which is similar to the dilog, but which only allows real solutions (the 
term Z7rln|i>fc| in Eq. (|8]) ensures this). In the chiral limit with M* > m + fh these two 
functions are equivalent. In addition they have the term (8 — 1/2) In |1 + 25\ whereas we 
have (5 + 1/2) ln|l + 25\. We believe this difference is due to a typographical mistake as 
the argument of the log function is linked to the multiplier outside, so that there should be 
no difference between them (the missing multiplication factor of two is easily accounted for, 
but not the sign change). 



A. Chiral Limit of P —* 77 

There is a well-known chiral limit of the P — > 77 case, namely n° — > 77 HH, when 
g^jy = e 2 Ncg n oQ 2 /4n 2 m. This implies that J = 1/2 for the pion. We can establish this 
from our full formulae. The chiral limit implies that M^o — ► so that the appropriate form 
of J is 

J = 2m 2 Li 2 (p,0)/M 2 , 

and we define M/m = e with e — > as M — > 0. We subsequently find p = e, cos0 = e/2 
and therefore 



2 1 rlnd-ex + x 2 ) 
e 2L 2 Jo x J 

-- f dx[-e + (1 - e 2 /2)x + e(e 2 + 3)x/3 + . . 
e A Jo 



= i-e 2 /12 + 6 4 /9 + ... 
Therefore J reproduces the 7r° — > 77 result in the chiral limit. 



VI. HEAVY QUARK EXPANSION 

Since we are particularly interested in the heavy meson decays D* — > D'j and B* — > S7 
we feel it is of interest to examine the heavy quark expansion of our loop integral J. To 
derive this we consider an expansion in terms of the light to heavy quark mass ratio in each 
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of the loop integrals. We make the arbitrary choice of m — m q and m = tuq, where m q is 
the light quark mass, and tuq the heavy quark mass. These lead to the following definitions: 

? = « (19) 
m 

M = m + A, (20) 

M* = rh + A* (21) 

where e — > in the heavy quark limit, and A, A* is the combined binding energy and light 
quark mass. 

A. Heavy Quark Expansion of J 

When dealing with J, Eq. (GO) and (0) lead to 



M/m = r + 1/e, 
M*/m = r* + 1/e 

where r = A/m and r* = A*/m which are independent of the heavy quark mass. We express 
J of Eq. (0) using these relations to find 

7 ft; I du [ dv- 



/o jo u/e 2 — 2(r* — r)uv/e — (2r + l/e)u(l — u)/e 

ignoring the constant term in e. This reduces to 

jHQL e(r* In |2r*e| — r In |2re|) 
(r* — r)(l + 2e(r* + r)) 

as the highest order terms in the expansion. Note that this term is of the form e In e and we 
can thus expect slow convergence of the heavy quark expansion. 

B. Heavy Quark Expansion of J 

We maintain our definition of e and r; however since J corresponds to J with m <-> m; 
we shall need the relations 
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from which 



M/fh = 1 + A/rh = 1 + re 
M*/m = 1 + A*/m = 1 + r*e, 



;i /■!-« l + (e-l)u 

J ps / cm / ai>- 



JO 



1 — u — 2(r* — r)ewt> — (1 + 2re)w(l — u) ' 



ignoring the e 2 contribution. Thus our heavy quark expansion of J reduces to 

JHQL = [Li 2 (l + 2r*e, 0) - Li 2 (l + 2re, 0) + 2e ^\~^ ( r * \ n |2r*e| - r In |2re|)l. 

To simplify this form further, we consider an expansion of the dilogarithm. Since 
Li 2 (l + 2r*e, 0) - Li 2 1 + 2re, 0) = - / 



,l+2r* e lull -d 
CtM 

r2re U 
r 2r*e 



rzr e 

- dz In z(l - z + z 2 - z 3 + . . .) 
-2e(r* In |2r*e| - r In |2re| - (r* - r) + 0(e)), 



then we arrive at our final form 



jHQL = j —( r * l n |2r*e| - r In |2re|). 

Once again observe the elne dependence, indicative of slow convergence. 

It appears that both expansions J H Q L and will only converge slowly to their true 

counterparts J and J. Thus, unfortunately, they are not so useful approximations for either 
the c or b quark cases. 

There are other possible expansions we could consider, namely that of the P — > 77 and 
V — > P7 loop integrals where there is only one quark flavor in the loop (such as r\ c — > 77, 
f]b — > 77, J ji> — > Vcl or T — > ^57) and consider some expansion as the quark mass becomes 
large. Unfortunately, such an expansion fails to be a good approximation, simply due to 
the assumption one has to make about the pseudoscalar and/or vector mass. For example, 
in the P — > 77 case one would assume the pseudoscalar mass M would consist of the sum 
of the quark masses along with some binding energy so that M = 2m + A. Then the loop 
integral could be expressed as 
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where p = M/m = 2 + A/m and one would attempt to do some sort of expansion near 
p = 2. Unfortunately such an expansion is impractical as the integral contains a pole at 
p = 2,u = 1/2. 



VII. RESULTS 

A. K* — > K^/ and the coupling ratio 

The observed K* branching fraction of 

T K*0—>K°'y/T K*+^K+-y = 2.31 i 0.29 

and corresponding coupling constant ratio of 

\g K ^ K 0j/9K*+K+y\ = 1-514 ± 0.125, (22) 

is far from its SU(3) predicted value of 2, but is simply understood in the quark loop 
formalism as shown by Bramon and Scadron ||. We quickly re-iterate this point. Using 



Eq. ([L5D and assuming g Vus = g Vds and g Pus = g Pds then 

g K *o KOl JdA K *°, K°}/m d + J s>d [K*°, K°}/m s 



-1.47, (23) 



g K *+K+ 1 2J U , S [K*+, K+]/m u - J S , U [K*+, K+]/m s 

where we have used a more complete notation, J q ,q>\V, P] to denote J for quarks q, q', vector 
meson V, and pseudoscalar meson P of masses m, rh, M* and M respectively. We used 
quark masses m u = m d = 340 MeV and m s = 510 MeV. The experimental uncertainty in 
the coupling ratio permits the s quark mass range: 475 < m s < 545 MeV, when m u = 340 
MeV, and a fixed s quark mass of m s = 510 MeV permits a u quark range of 225 < m u < 385 
MeV. Since an s quark mass of m s = gives such a good comparison between the quark 

triangle diagram and the experimental measurement, we will choose such a mass throughout 
this work, along with m u = m d = 340 MeV. The result of Eq. (|23| ) compares well with that of 
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|| , indicating their chiral limit formulae is appropriate. In fact we can observe the variation 
from the chiral limit to M — Mk using our J; as Figure |2| shows there is very little change 
in the result. 

We can also dramatically show how the s quark mass breaks the SU(3) symmetry. Fig- 
ure |3| displays the behaviour of gK*°K y/ 9k*+k+'j from m s = m u = m d = 340 MeV (the 
SU(3) limit) to m s = 550 MeV clearly indicating that it is the violation of constituent 
quark masses from SU(3) symmetry that is responsible for the large deviation in K* mesons 
from the expected symmetry. The K* radiative decays are particularly sensitive to SU(3) 
violations as they involve the constituent masses of strange and non-strange quarks in the 
loop of the corresponding triangle diagram. Subsequently, the heavier meson cases are best 
computed by the triangle diagram which can allow for different constituent masses in the 
loop, rather than an SU(4) or SU(5) symmetry. 



B. Measurements of 



9Pqq' 



We may obtain estimates for the gp q q/ coupling constants using experimental measure- 
ments of the P — > 77 decay widths. In particular the widths for tt° — > 77, 77 — > 77 and 
r/f — > 77 processes are well known and we should be able to determine gp ua and gp sS to 
reasonable accuracy (we make the isosymmetric approximation gp ua = gpddi m u — md)- 
Beginning with the 7r° meson which is the antisymmetric mixture (uu — dd) / v^2, we find 

6e 2 {/2\ 2 J u [n°] 



<?7T°77 /l„2 1 Q / 9-k°uu 



e 2 \ a 9p^_9p^\j ^ i 







(-1 




1 m d J 



6vr 2 1 ^2 y/2 



e 



2 



gpuuJu[^°}/m u . (24) 



2V27T 2 

In a similar fashion, we use the standard octet-singlet pseudoscalar mixing angle 9p to 
ascribe the 77-77/ mixing by 

77 = — ^[(cos^p — V2sin6p)(uu + dd) — v^sin^p + cos 6 p)ss], 
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r\l = —i=[(sm6p + \/2cos9p)(uu + dd) + >/2(cos Q P — y/2 sin 9p)ss\. 
V6 

Following a methodology like Eq. (E3j) we obtain relations between the covariant amplitudes 
and meson-quark couplings which are given in Table |. We have used m u = 340 MeV, 
m s = 510 MeV, two different mixing angles, Op = —10.5° (in accordance with the quadratic 
Gell-Mann-Okubo relation) and 9p = —20° (which is partly favored by pseudoscalar decay 
processes) along with the covariant couplings from the measured decay rates |7j 

T P ^ 77 = m P g 2 P ^/64ir, 

to obtain several estimates of the pseudoscalar-quark couplings as given in Table ||. We 
used the rj — > 77 and rj/ — > 77 to determine simultaneously the values. 

From the results, it appears gp U u differs as determined from ir°, rj and rj/ processes. Con- 
sidering that the 77 meson is about four times as massive as the pion, it may be appropriate 
to allow for such a mass dependency in the coupling constant. Suppose we label the first 
coupling constant from 7r° — > 77 as gp U u (rn^ ) , while the second from 77 — > 77 and rjf —>■ 77 
as a coupling constant somewhere between m v and m n ,. Numerically we took the appro- 
priate mass as the equal weight average, (mi + mL)/2. By linearly interpolating between 
these two couplings, we estimated a value of gp U u( m ^) = 4.61 ± 0.19 for Op = —10.5° and 
gpuu(m 2 v ) = 4.39 ± 0.17 for P = -20° . 

The Goldberger-Treiman (GT) relation at the quark-level, gives us a good check of our 
results. For the pion, the relation reads 

U9Puu(mlo)/V2 = m u . 

Using our coupling value in Table along with m u = 340 MeV we predict f n = 93.5 ± 3.46 
MeV which compares well with the experimental result f n = 92 A ± 0.26 MeV [0. 

Also included in Table [TJ is the estimate of gp c5 using a charm quark mass of m c = 1550 
MeV along with the experimentally determined width [0] of r Vc ^, 71 = 7.0 ± 2.6 keV. 
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C. Measurements of gv q q' 



There exist many useful decay channels V — > P7 and corresponding data from which 
we can determine the product gvqq'Qpqq'- To this end we proceed in two steps. Firstly, 
we interpret individual meson-meson-photon couplings in terms of meson-quark-antiquark 
couplings, deriving relations between them as shown in Table |. Assuming isospin symmetry 
there are only two unknown products of couplings involved in the light meson sector; one is 
9vuu9puu for non-strange quarks while the other is gvssQPss for strange quarks. 

Following this we extract individual meson-meson-photon couplings from the most re- 
cently measured decay widths Tv-,p-y = {ni v — mpYg^ pry /12TTrriy by simply removing the 
kinematic factors. The results are listed in the first column of Table |TJ. As one can see, they 
scatter over a relatively wide range. 

We are able to determine gvuudPuu solely from any one of the processes p° — ► 7r°7, p + — > 
7r + 7, p° — > 777, uj — > 7r°7 and uj — > 777. In addition, the decays uj — > 777 and — > 777 can 
be used to simultaneously solve for gvuugpuu and gvssgpss- Our numerical results are shown 
in the second column of Table [H| where we use the same quark masses as previously along 
with standard mixing angles. 

The values of the product gvuugpuu turn out to lie in a quite small range, except for that 
from the — > 7r°7. However, it would fall into this range had we chosen a mixing angle 
of about By = 224°, a change of 4.6°. Such a high sensitivity of <p — > 7r°7 to change in 
mixing angle suggests it is reasonable to exclude this channel from our analysis. Recalling 
the couplings of pseudoscalar meson with quark-antiquark pairs discussed previously, we 
now obtain gv U u- 

As the light vector meson masses vary by less than 30%, we shall not attempt to distin- 
guish between the slightly different couplings and thus on average we find gvm — 2.40 ±0.08 
(weighted average) for 6p = —10.5° and gy U u = 2.35±0.08 (weighted average) for 6p = —20°. 
It differs from gp u u, revealing a substantial violation of the spin symmetry in the triangle 
scheme. 
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We repeat this procedure in the analysis of gvss, but with fewer channels to determine a 
result. Consequently we have gv s s = 1-10 for 9 P = —10.5° and gv S s = 1-06 for 6 P = —20°, 
indicating a large SU(3)y symmetry breaking once again. Estimates for gv C c using the 
J/ijj — > ^ c 7 channel yield gv C c = 0.92 ± 0.23. Note that gy C c and gp c5 are not substantially 
different, perhaps indicative of a limit gy qq - = gp g g as m q gets large. 

For completeness, we wish to obtain a measure of gv u s using the product gvusgpus- 
However, we have no means of getting gpds for the kaon in the triangle scheme. This is 
because unlike 7r°, the K° — > 77 decay is not mediated by pure electromagnetic interactions. 
However assuming the Goldberger-Treiman relation at the quark level 

fKgpds(fn 2 K ) = (m u + m s )/2 

we find gpds = 3.77 ± 0.03 where we have used fx = 113.0 ± 1.0 MeV 0. Subsequently, 
gvds = 2.21 ± 0.10 (averaged over the charged and neutral processes). 

We ought to point out that the triangle scheme has also been applied to radiative decays 
of 77/ into p° or u, but failed to yield coupling constants near the above range. This suggests 
to us that we should treat r]/ in a different way which would most likely incorporate the U(l) 
anomaly. 

VIII. PREDICTIONS 

A. — > 77/7 coupling constant and branching fraction 

We can use our best fit estimates of the meson-quark coupling constants to predict the 
decay width for the decay (f> — > r/fj. 

g^ nh = -^{(cos6 v - v / 2sin6 l y)(sin6 l p + V2 cos, 9 P )g Vm g Pm J uA^^ V']/ m u + 
2(sin#y + v / 2cos6'y)(cos^p - \/2sm6p)gvssgpssJsA < i ) i V']/ m s} 
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and using m u = 340 MeV, m s = 510 MeV and mixing angles 6p = —10.5°, 9y = 219.4° along 
with our couplings from Table we compute the coupling to be g^- vr p 1 = —5.69 x 10~ 4 . 
This gives a branching ratio of 

Br((f) -> 77/7) = 4.16 x 10~ 4 

which is slightly above the experimental upper limit of Br(<f) — ► 77/7) < 4.1 x 10~ 4 at 90% 
confidence level 0. However, we note that a change to an s quark mass of m s = 500 MeV 
produces a branching fraction of Br(<f) — > 77/7) = 3.24 x 10~ 4 , so that the result displays 
very sensitive dependence on the choice of s quark mass, and probably vector mixing angle. 
Also we remain cautious of predictions involving the 77/ meson due to its association with 
the U(l) anomaly. 



B. D* — > D7 and B* — > .67 coupling ratios 

Since our J is approximation free (i.e. no chiral limit assumptions) we can safely use it 
in the D* and B* meson cases. We do assume gy u Q = g Vd Q and gp u Q = gpdQ where Q is 
either the c or b quark (much like we did in the K* — > K'-f case) to obtain 
g D *o DOl 2(J U , C [D*°, D°}/m u + J c , u [£>*°, D°]/m c ) 



and 



(?d*+d+ 7 - Jd, c [^* + , D+]/m d + 2J Cid [L>*+, D+]/m c 



(25) 



(26) 



<?b*+b+ 7 ^+]/m u - 2J b:U [B*+, B+]/m b 

Relations (^) and ( p6| ) allow us examine the coupling constant ratios as a function of the 
c and b quark mass, respectively. The results appear in Figures |] and [|. In order to 
give actual values we use a c quark mass of 1550 MeV (approximately half the J/ip mass) 
yielding gD*°D -y/gD*+D+- r = 6.47 and a b quark mass of 4730 MeV (approximately half the 
T mass) which gives gB* B°-y/ gB*+B+-y = 0.018. We can compare our results with those of 
other workers. These are presented in Table JTJ. We hoped that our study of gvqq',gpqq' 
measurements would enable us to make some reasonable guesses of gvucgpuc and gv u l9pubi 



but the data does not allow this. Thus we cannot make predictions about actual decay 
widths. 



IX. CONCLUSIONS 

We have successfully evaluated V — > P7 and P — > 77 processes in a quark triangle 
diagram scheme which is valid for arbitrary vector or pseudoscalar masses. By comparison 
with available experimental data, we found that this scheme works well for all radiative 
processes involving the light mesons (no charm or bottom quarks), except for — > 7r°7 (due 
to the sensitivity of this channel to the mixing angle) and 77/ — > p°(cu)^. 

The scheme produces well determined estimates of the meson-quark-antiquark couplings 
for the light mesons. The large difference between g Vq q> and gp qq -i indicates a substantial 
violation of spin symmetry in the quark triangle formalism. We also observed a relatively 
weak SU(3) chiral symmetry breaking due to the finite masses of the Goldstone-type pseu- 
doscalar mesons, along with a more apparent SU(3)y symmetry breakdown arising from 
the difference in light constituent quark masses. We note that these coupling constants are 
relatively insensitive to the pseudoscalar mixing angle. 

A number of predictions have been made based on the scheme. Firstly we note that 
our theoretical result for the — > 77/7 decay width is around the present experimental 
upper limit and awaits comparison with further measurement. Secondly our prediction for 
5 , D*°D°7/fi , D*±z)±7 = 6.47 with m c ps mj/^/2 is within range of other theoretical estimates, 
while gB*°B°y/ 9B* ± B ± -y — 0.018 for mb mx/2 is small compared with the few results in 
the literature. We expect future measurements of these radiative decays will distinguish 
between these predictions. 
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Process 



TABLES 

Relation between covariant couplings and meson-quark couplings 



7T U 

V 



Vc 




p 
p 
p' 

CO 







77/. 

7)1 







K*° - 
K*+ - 



► 77 
77 



77/ — > 77 



' 77 

7T°7 



7T + 7 



?77 



7T°7 



A 



^7 



7T°7 



-> 777 



,2 



_ e" 



m 



9vrr = ^6^{ 5 ( cosdp ~ V% sin6 P )g Puil J u [ri}/ m u - 
V2(sin9 P + V2 cos6p)g Ps sJ s [i]]/m s } 
■{5(sin 6>p + V2 cos 9 P )g Pua J u [r]f}/m u + 



^ V S - , „- -L.J, 

>/2(cos P - v 7 ^ sin 9p)g Ps s J s M /m s } 

.2 _ r . . 



,0_ >7r 7 — 4^2"5'V'MMfl , PMM^M,u[p ! 71 " ]/' m u 



dp+^ir+y = -^2 9vud9Pud-Ju,d[p + ,K + ]/m u 
g p o^ rr/ = ^§(cos9p - sin 9 p)g Vu ugpuuJu,u[P iV]/ m u 
g^o-y = ^J(sin6» y + \^cos9 v )g Vu ugpuuJu,u[ UJ ,n°]/m u 

9w->rn = 12^2 { 

(sin9 v + v / 2cos6'y)(cos 9 P - \/2 sin 9 P )g VuU g PuU J u ^ u {uj, r]\ /m u + 
2(cos6» y - \pl sin 9 V ) (sin 9 P + N/2cos0p)#y S s#p S sJ S)S [u;, ?7]/m s } 
Sty^o 7 = iJ( sin6l p + cos P ) g v uugpuu J u ,u [v', P°]/ m « 

= 19^2 { 



9r)t->ury 1271 

(sin#y + cos (9y) (sin 6> P + cos 9p)gvuu9PuuJu,u[v f , u}/m u - 
2(cos(9 y - v 7 ^ sin 9 V ) (cos 6> P - ^/2sin9 P )g VsS g PsS J StS [r]f, u]/m s } 
g^^ = (cos 9 V - V2sm9 v )g VuU g Pu uJu,u[ ( l ) ^ ]/ m u 



~ 12n 2 I 



iT°7 
K+7 

► ^7 



(cos 6V - v 7 ^ sin 9 V ) (cos 6> P - V2 sin 9 p)g VuU g PuU J u>u [(p,r]]/m u - 
2 (sin 6V + \[2 cos 9 V ) (sin 6» P + v / 2cos6^)#y s jgp sg J SiS [0, ^/V,} 



g K *+->K+~i = -£?gvu7sgp U 7s(?J UyS [K* + , K + ]/m u - J S , U [K* + , K + ]/m s ) 

9J/ip-*vcj 



= -^9vcc9PccJc,c[ J /^^Vc]/m c 



TABLE I. Radiative decays of ground state mesons and relations between covariant couplings 
gp 77 ,<yVp 7 or gpv-y and meson-quark-antiquark couplings gvqq>,9Pqq> 
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Experimental result Meson-quark-antiquark coupling 11 

:10- 4 MeV _1 ) Q P = -10.5° 6 P = -20° 



|^o 77 | = 0.2516 ±0.0091 
\g mi \ = 0.239 ±0.011 



''77 I 



0.312 ±0.016 



gpuu 

9Pss '- 



g Pu u = 5.U±0.19 
: 4.03 ±0.14, J g PuU = 3.56 ±0.12, 

6.42 ±0.67 g PsS = 8.19 ±0.65 



I5V77I 


= 0.07297 ±0.01366 


g Pc5 = 2.03 ± 0.38 


Iv- 


^o 7 | = 2.96 ±0.38 


gvuugpuu = 14.82 ± 1.90 


l#p+- 


^+ 7 | = 2.24 ±0.13 


9vuu9puu = 11-32 ± 0.64 


Iv- 


^ 7 | = 5.67 ±0.53 


gvuugpuu = 10.96 ± 1.02 gvuugpuu 




>7r o 7 | = 7.04 ±0.21 


gvuugpuu = 12.53 ± 0.38 



Ift^l = 1.83 ±0.23 
|^^ 7 | = 2.117 ±0.052 
|^^o 7 | = 0.417 ±0.021 
|#k*o^o 7 | = 3.84 ±0.17 
|fe*+^+ 7 | = 2.534 ±0.115 



9Vuu9Pu 



12.3 ± 1.5, 



9Vuu9 Puu 

gvssgpss = 7.08 ± 0.17 I gvssgpss ~- 
gvuugpuu = 25.1 ± l.3 b 

9vds9Pds = 8.43 ± 0.37 



= 9.56 ±0.89 

= 10.7 ± 1.3, 
8.66 ±0.21 



\9j/iI>- 



1.67 ±0.26 



9vus9P-b 



9Vcc9Pc 



8.21 ±0.37 



1.87 ±0.30 



TABLE II. Determination of meson-quark-antiquark couplings 
m u = m d = 340 MeV, m s = 510 MeV, m c = 1550 MeV, V = 219.4° 



gvuuQp-L 



11.9 ± 0.6 for By = 224° 
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Reference ^d^d^^I 9d*+d+~,\ \9b*°b°j/9b*+b+j 

i 

This paper 6.47 0.018 

@ 3.05 ± 0.63 0.49 ± 0.38 

[0] 2.98 ± 0.62 

@ 6.32 ± 2.97 0.64 ±0.51 

||] a 11.0 ±2.3 

[|l| b 12.9 ±2.7 

3.28 ± 0.83 
@ 5.4 ±1.1 

- 0.58 ± 0.44 

@ c 6.95 ±1.44 

@ 6.61 ± 1.37 0.62 ±0.48 

@ 5.54 ±3.00 0.59 ±0.48 

[0 60 0.8 d 

@ 3.50 ± 0.73 

@ - 0.678 ± 0.523 

H 3.84 ± 0.80 

H c 1.66 ± 0.34 

|3 f 3.93 ± 0.84 

|U s 4.49 ± 0.96 

[|| h 3.92 ± 0.84 

H 2.00 ± 0.41 

[|| 3.78 ± 0.78 

TABLE III. Summary of theoretical estimates 

a Gaussian wave function with m u =300MeV, m c =1500MeV 
b BS wave function with m u =350MeV, m c =1500MeV 
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c m u = m d = 300 MeV, m s = 500 MeV 

d The original paper contained an error for this calculation. This is the corrected value. 

e Zero anomalous magnetic moment of charm quark 

f SU(4) symmetry 

^Broken SU(4) by Ml transition 

h Broken SU(4) by 1/Af£ 
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FIGURES 




FIG. 1. Quark triangle diagrams contributing to V — > P7. 
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9 K *+K+ 1 




m K (MeV) 



FIG. 2. Appropriateness of chiral limit, shown in the lack of sensitivity of 9K*°K°~f/ 9k*+ K+-y 
to K meson mass. 
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m s (MeV) 



FIG. 3. Breaking of SU(3) by s quark mass. The experimental measurement is included. 
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g D *o D o 



9d*+D+j 
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